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We analytically establish the role of a spectrum of Lyapunov exponents in the evolution of 
phase-space distributions p(p,q). Of particular interest is A2, an exponent which quantifies the 
rate at which chaotically evolving distributions acquire structure at increasingly smaller scales and 
which is generally larger than the maximal Lyapunov exponent A for trajectories. The approach is 
trajectory-independent and is therefore applicable to both classical and quantum mechanics. In the 
latter case we show that the h — > limit yields the classical, fully chaotic, result for the quantum 
cat map. 
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The canonical measure of chaotic Hamiltonian dynamics is the non-zero Lyapunov exponent jjj which quantifies 
the rate of exponential divergence of neighboring classical trajectories. This measure has resisted translation to the 
. Hilbert space of phase space distributions, where the linearity of the evolution equation precludes the asymptotic 
exponential divergence of initially "close" distributions. Rather, chaos in distributions is described in terms of the 
character of the spectrum of the Liouville operator Q), a difficult property to access. As a consequence, a quantitative 
I j diagnostic for the manifestation and characterization of chaos in the evolution of distributions has proven elusive, 
. although several conjectures and limited diagnostics have been proposed 

In this article we analytically establish a diagnostic of chaos for distributions that is independent of trajectories 
and is valid for arbitrary initial distributions and at asymptotic times. In particular, we expose the role of an entire 
spectrum of generalized Lyapunov exponents A 7 || in the evolution of distributions. Although we derive the results 
by considering phase-space trajectories, our final expressions for A 7 depend solely on the evolution of the phase space 
j distribution. We thus have a means of studying the role of Lyapunov exponents in the evolution of both classical and 
ON ■ quantum systems, as shown below. 

We also show that theparticular exponent A2, which is generally larger than the standard trajectory-based maximum 
Lyapunov exponent A |8[, assumes particular importance as the arbiter of the growth of phase space structure. This 
growth has been thought fl to characterize chaos for distributions. In particular, we show that (1) a distribution 
evolves so as to acquire structure at increasingly smaller scales at an exponential rate given by A2 (2) the rate at 
which the information in the system moves to the smaller scales (and hence, the rate at which it relaxes) is given by 
. A2 and (3) the error in a calculation for a fixed level of resolution increases exponentially in time with A2. 

Consider the usual computation of the Lyapunov exponents |lj : Let the equations of motion of a point in phase 
space be Xi — fi{x), where x[= (q,p) for Hamiltonian systems] is the vector denoting all phase space variables and 
individual vector components are denoted by subscripts. Our variables are scaled to be dimensionless so that we 
may define a metric in phase space. The equations of motion for the vectors in the tangent space are obtained by 
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linearizing the equations of motion [substituting x + q into xi = fi{x)] around this fiducial trajectory: 

j 

where My = ^j 2 - is the Jacobian matrix at the point x(t). If this system is evolved forward in time, then we can 
obtain the maximal Lyapunov exponent as [Theorems A and B, Ref. jij] 

limi]n(k(*)l) = M*) (2) 

t— »oo t 

where • • • | defines the standard norm in phase-space and the result is independent of <; except for a set of vanishing 
measure. For a Hamiltonian system the set of exponents is unchanged under the transformation A — » —A as a 
consequence of time-reflection symmetry. 

We use the foregoing to derive the equations of evolution for the first-order derivatives of a distribution. The 
assumptions necessary in our analysis imply that the results are valid for all distributions that evolve non-trivially 
under the Liouville equation and are differentiable: This specifically excludes delta-function distributions (point 
trajectories) and time-invariant solutions to Liouville's equation. Consider the relationship between the densities at 
the point x and x + q: 
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p(x + <r) = p{x) + ^-q^p( x ) 



(3) 



where we neglect higher-order terms, as in the trajectory-based derivation. Taking the total time-derivative of both 



sides and noting that 
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where here and below p and its derivatives are evaluated at x. Using Eq. ([!]), we get 
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the last equality following from the independence of the various • Eq- (Q) and Eq. (|7|) for the evolution of £ and Vp 
respectively are identical except for a minus sign. This makes physical sense: Since the number of points cannot be 
created or destroyed along trajectories, the density sharpens (flattens) along the direction in which points move closer 
(further). This being the case, A is directly reflected in the gradient of the distribution Vp(x(t)) at asymptotic times 
at almost every phase space point JiTj[ | in the chaotic region. Remembering the reflection symmetry of the Lyapunov 
exponents, the maximal Lyapunov exponent may hence be equally well computed as 



X(x)= lim -ln(|Vp(*(t))|). 
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We now consider the averaged quantites x-y where 
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(8) 



(9) 



Here Tr indicates the trace or integration over phase space and 7 is an arbitrary real number; note that the denominator 
Tr[p 7 (x(£))] is a constant for Hamiltonian evolution. Several interesting properties have been established for similarly 
constructed averages using ? instead of Vp 0. Using the same approach, since Eqs. (||) and (JsJ) are identical, we get 
that the following properties also hold for A 7 defined via Eq. (||): (1) A spectrum of generalized Lyapunov exponents 
can be defined through 



lim -ln(x 7 ) 
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(2) The 'usual' maximal Lyapunov exponent A is a member of the spectrum: lim 7 ^o A 7 = A and (3) A 7 is concave in 
7; i.e., if 71 < 72, then A 7l < A 72 . Note that by defining the averages in terms of Vp rather than the quantities x 7 
and A 7 are independent of trajectories. 

We note that the dependence of A 7 on 7 arises due to the variation in local stretching rates in different regions of 
phase space. Hence, for a system where the local stretching rate is a constant, A 7 is independent of 7 and equals the 
maximum Lyapunov exponent A. 

Considerable physical insight emerges by considering A2 and the quantity X^it) i n Fourier space. That is, consider 
the Fourier expansion of a distribution p(q,p) — m Pn,m(t) fn,m where f n ,m{p,q) = cxp{27ri(rtp + mq)}. We use 
a two-dimensional discrete Fourier basis for simplicity; all arguments generalize to multiple dimensions and Fourier 
integrals. In Fourier space X2(t) is of the form 



_ Tr[ \Vp(x(t))\ 2 ' 
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Hence xi is t ne root-mean-square Fourier radius, measuring the Fourier space extent of the phase space distribution. 
Equation ( |i"o| ) shows that 
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lim iln( X2 )=A 2 ; (12) 
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or that, for a chaotic system, \2 increases exponentially with time. Thus, since the higher (\n\, |m|) modes correspond 
to structure at smaller scales, a distribution moves exponentially in time, with rate A2, from structure at a discernible 
scale to structure at extremely small scales. We can relate this to the loss of accuracy associated with the chaotic 
evolution of point trajectories as follows Jll| . All the information about the distribution is encoded in the initial Fourier 
basis expansion; as \2 grows, this information leaves any finite range of |n|, \m\ exponentially fast. Attempting to 
account for all Fourier modes corresponds to retaining an infinite amount of information and is inconsistent with a 
finite-resolution measurement in phase-space. Thus, Fourier modes with mode-numbers greater than some n max , m max , 
where — — , ^ — are the limits of resolution in p, q respectively, are effectively non-observable and the information 
therein not retrievable. In this sense the evolution of a chaotic distribution can be reconciled with the exponential 
loss of accuracy inherent in the chaotic evolution of trajectories, also a consequence of finite measurement capability. 

Equation ( [To| ) provides deep insights into the relation between the Lyapunov exponents and distribution dynamics 
but should not be regarded as a competitive tool to compute Lyapunov exponents when trajectories can be calculated. 
However, since x~t(t) 110 longer relies on trajectories, we may now study it for quantum mechanical phase-space 
distributions, thus analyzing "quantum chaos", i.e. the effect of chaos on quantal evolution. Again, \2 proves 
particularly valuable since, by integration by parts 

a _ Tr[ [Vpj 2 ] _ Tr[pV 2 p] 
X2{ > ~ 4vr 2 Tr[p 2 ] " 4^ 2 Tr[ (0 2 ]- [ ' 

Noting that V 2 p = {p, {p, p}} + {q, {q, p}}, where { , } denotes the Poisson bracket, we may now perform the standard 
quantization {A, B} — > — i[A, B] to yield a quantal \i which is independent of the representation. In the Wigner-Weyl 
p3[ representation it has precisely the classical form [Eq. (|l3|)], with p replaced by the Wigner function W(q,p). 

As an example of the insights afforded by this approach we consider the cat map on the torus, which is a classical 
K-system ]l2] ] and which has been recently shown |Q to display smooth quantum-classical correspondence. The 
dynamics of this system, whose classical and quantal propagators for distributions are known analytically p4| , are 
those of a kicked oscillator with Hamiltonian [[flj 
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H=p 2 /2p + eq 2 /2 6(8 -t/T). (14) 
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restricted to a torus < q < a, < p < b. The cat map corresponds to the choice r\ = Tb/pa — 1 and £ = —eTa/b = 1 
and a — h/ab acts as a dimensionless form of Planck's constant for this problem. Two other cases of interest we 
study are the (stable) elliptic map (ry = 1,£ = —1) and the borderline parabolic case (77 = 1, £ = 0). We note that 
this classical system has a constant Jacobian matrix and x-y = X is independent of 7. 

We then anticipate that (1) for the stable system x oscillates as a function of time, (2) x grows linearly for the 
parabolic case and (3) a chaotic system has x growing exponentially in time with the asymptotic rate given by A. 
This behavior is indeed confirmed computationally, as shown in Fig. 1, for these three maps. We emphasize that (1) 
the behavior is independent of the initial distribution; the distributions used here were randomly initialized and (2) 
in the chaotic case, there is no saturation of x(t) even though the phase-space is bounded. 

In the quantal case we study the dependence of X2 (t) on the degree of classicality of the system, that is the value of 
a. To eliminate kinematical effects of changing a we consider the evolution of Wigner distributions which are initially 
of the Gaussian form 



W(q,p) = Npexp[ _ — ]exp[ _ — ] (15) 



where (3 is a scaling variable (see below), Np is a normalization factor, (/3qo, /3po) specifies the location of the Gaussian 
of width (/3<7g, f3<J p ) in a phase space of dimension f3 2 ab. Since a — h/ f3 2 ab, by increasing (3 we approach the classical 
limit while preserving the ratio (/3 2 a q a p / /3 2 ab) of the volume of the initial distribution to the volume of phase space. 
Here we choose this ratio as w 0.0175, sufficiently large to display quantum effects at relatively early times. Note 
that for these large initial distributions, the approximate methods of measuring the separation of centroids or the 
growth of second moments fail to detect the existence of chaos, even for the classical system. Our results are 
displayed in Fig. 2. 

For the classical cat (h = 0,/3a = f3b = 1), shown again for comparison, X2(t) grows exponentially after a short 
transient, with a rate equal to A. The behavior of %2 for the quantized cat map is shown for a — 10~ 5 , 10~ 2 and 
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10 _1 . Near the classical limit (a = 10~ 5 ;h = l,(3a = j3b= 316.2) \2 exhibits exponential growth, indistinguishable 
from the classical behavior on the finite Fourier grid (of size 2048 x 2048) used in our computations. However, as a 
is increased to 10~ 2 (h = 1, (3a = (3b = 10) the initial growth of \2 is seen to be faster than that of the classical cat 
map, suggesting that the amount of structure in the quantal distribution at short times actually increases relative to 
the classical one. This is indeed the case since one expects the quantal W(q,p) to adiabatically follow the classical 
p(q,p) but with added fringes |u| (interference structures). Such interference structures can be seen in contour maps 
of the evolving distributions. For example, an examination of Figs. 1-3 of Ref. JfJ at t — 4T shows clearly that the 
distribution with a = 10 -2 has far more structure at finer scales than does the distribution with a = 10~ 5 . This, 
of course, implies greater support for the distribution at larger Fourier node-numbers, leading to larger \2 values. 
However, supra-classical growth does not persist for larger a or for longer time because a second quantal effect kicks 
in: Quantal distributions resist the growth of structure at scales smaller than fi |l5[ |, implying that the support of 
the distribution travels slowly across the — boundary in Fourier space. This is marked by a clear slow-down of the 
growth at % 2 ~ 100 = 0{— ). As a is increased to 10 _1 (h = 1, f3a = (3b = 3.162), the initial rapid growth saturates 
in one time step, resulting in a slowly growing \i- The interplay between the two effects thus implies in general a 
finite value of a at which the quantum system acquires structure maximally rapidly for early times. Note also that 
the initial differences between the various values of \i m Fig- 2 is exaggerated by the fact that we are using discrete 
time-steps which are themselves of the order of in our model system. 

Since the quantal propagator is a smooth function of a, the results at larger values of a imply that on the full 
(infinite) Fourier grid the quantal \2, even at the smallest values of a, will ultimately slow down, reflecting the 
non-ergodic nature of the quantum map |14|. However, and this point is significant, for any finite resolution there is 
always an a sufficiently small such that the classical and quantal behavior are indistinguishably chaotic. Thus, the 
various govern the behavior of classical and quantal distributions in precisely the same manner in the limit h — > 0, 
as manifest in the behavior of X2- Note that the form of x%(t) enables us to verify this behavior without masking 
effects arising from the choice of broad initial distributions or from saturation and relaxation in the finite phase-space 
(see the discussion in Ref. JtJ). 

In summary, we have provided a quantitative demonstration of the role of Lyapunov exponents in the behavior of 
classical and quantal distributions and have identified \2 (t) as a measure of the increasingly fine detailed structure of 
the distribution. We have thus related the generalized Lyapunov exponents directly to the properties of distributions, 
independent of trajectories. Our approach is valid for arbitrary initial distributions and asymptotic times, unlike other 
methods ||[7|], and has been explictly derived from the definition of Lyapunov exponents, unlike previous conjectures 
In particular, it avoids the need to introduce perturbations to define the stability of distributions Q and 
avoids basis-dependent definitions of chaos in distributions 0. Rather, it provides a consistent understanding of the 
role of Lyapunov exponents in the dynamics of classical distributions and affords a method of connecting Lyapunov 
exponents to the behavior of quantal distributions in the semiclassical limit % — ► 0. 
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FIG. 1. X as a function of time for the classical elliptic, parabolic and cat maps. 



FIG. 2. Xz(t) f° r the classical (a = 0) and quantal (a = 10 1 ,10 2 and 10 5 ) cat map. 
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